具体内容及详细教程请关注微信公众号:夫也的笔记
最近老师在课上安排了一个汇报作业,要求对某一地区的多年NDVI数据进行趋势分析,对于小白的我来说,真是!唉!但是功夫不负有心人,在我充分利用百度搜索、谷歌搜索和Bing搜索后,终于让我找到了相关教程和知识。只是需要用Matlab软件进行一系列操作,废话不多说上代码。
趋势分析有很多种,今天介绍的是Theil-Sen median趋势分析,一般还要结合Mann-Kendall建议,这个MK检验后续文章会写。
Theil-Sen median趋势分析是一种稳健的非参数计算方法,计算式如下
公式中SET指计算n(n-1)/2 个数据组合的斜率的中位数,ETi和ETj代表i和j年的ET值。如果SET>0,则ET呈上升趋势,否则,ET呈下降趋势。
Matlab具体代码如下:
% @author 1154318421@qq.com
[a,R]=geotiffread('E:\GIS\NDVI\2000.tif');%先导入投影信息
info=geotiffinfo('E:\GIS\NDVI\2000.tif');
[m,n]=size(a